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ABSTRACT 

iPatat et aTJ recently inferred the existence of circumstellar material around a normal Type la supernova (SN) for the first time, 
finding time-variable Na I D absorption lines in the spectrum of SN 2006X. We present high-resolution spectroscopy of the bright 
SN la 2007 af at three epochs and search for variability in any of the Na D absorption components. Over the time range from 4 d 
before to 24 d after maximum light, we find that the host-galaxy Na D lines appear to be of interstellar rather than circumstellar 
origin and do not vary down to the level of 18 mA (column density of 2 x 10 11 cm -2 ). We limit any circumstellar absorption 
lines to be weaker than ~ 10 mA (6 x 10 10 cm" 2 ). For the case of material distributed in spherically symmetric shells of radius 
~ 10 16 cm surrounding the progenitor system, we place an upper limit on the shell mass of ~ (3 x 10~ 8 )/X Mq, where X is the 
Na ionization fraction. We also show that SN 2007 af is a photometrically and spectroscopically normal SN la. Assuming that the 
variable Na D lines in SN 2006X came from circumstellar matter, we therefore conclude that either there is a preferred geometry 
for the detection of variable absorption components in Type la supernovae, or SN 2007 af and SN 2006X had different types of 
progenitor systems. 



Subject headings: circumstellar matter 
2007af) 



supernovae: general — supernovae: individual (SN 2006X, SN 



1. INTRODUCTION 

Type la supernovae (SNe la) are currently the only dis- 
tance indicator that can be used effectively out to cos- 
mological distances (e.g.. iRiess et alJI 1998t IPerlmutter et alj 
[UMlLeibundgut 2004; FiliD penkoll2005F Ries s et al .1 120071) . 
Therefore, understanding the nature of these explosions 
and any potential systematics that may be present is of 
great importance to cosmology. However, we still do 
not know what the progenitor systems of SNe la are, 
and observations suggest that there may be at least two 
physi c ally different progenitor classes (e.g., Mannuc ci et alj 
2005t IScannapieco & Bildstenl 12001 iMannucci et all 120061 
Sullivan et al.l l2006: Ouimbv et al. 2007), as well as some pe 
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culiar object s (e.g.. iLi et alJl200U 120031: lHamuv et al.ll2003l 
al though s ee Benett i et alJl2006h . 

IPatat et aTJ (l2007al) have recently made a possible break- 
through in the study of SN la progenitors by optically detect- 
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ing circumstellar material (CSM) in a SN la for the first time. 
Using high-resolution spectra of SN 2006X spanning from 
just before maximum light to four months later, they showed 
that at least four distinct components of the Na I D absorp- 
tion lines varied with time until ~ 2 months post-explosion. 
Although similar behavior has been seen in Milky Way stars 
and is generally attributed to smal l interstellar clouds moy - 
ing across the line of sight (e.g. JWertv & Fitzpatrickl 2001). 
in SN 2006X the lack of time evolution in the corresponding 
Ca II H & K absorption features at the s ame veloci ties proba- 
bly rules out that interpretation. Instead. IPatat et all conclude 
that the variable absorption is from circumstellar clouds in 
the progenitor system that were ionized by the radiation from 
the supernova and recombined several weeks later; because 
Na I has a much lower ionization potential than Ca II, the 
Na I line profiles can change without an accompanying effect 
in the Ca II lines if the ionizing radiation has an appropriate 
spectrum. These results appear to indicate a single-degenerate 
progenitor for SN 2006X with a red-giant companion. 

Multiple-epoch high-resolution spectroscopy is available 
for only one pre vious SN la, the peculiar SN 2000cx 
dPatat et alj|2007bl) . so it is not yet known whether the time 
evolution seen in SN 2006X is common. Assuming that the 
variable absorption is related to material from the SN progen- 
itor, if the behavior of SN 2006X is not universal then either 
there must be geometric effects that limit the visibility of the 
absorption to certain lines of sight (e.g., near the orbital plane 
of the progenitor system) or there are multiple progenitor sys- 
tems for SNe la. 

In this Letter, we present high-resolution spectra of the 
Type la SN 2007 af obtained at -4.3, +16.6, and + 23.7 d rela- 
tive to m aximum light. We use these data to test the lPatat et alJ 
(2007a) model, searching for variability in the Na D absorp- 
tion features. 

2. OBSERVATIONS AND DATA REDUCTION 

SN 2007af was discovered by K. Itagaki on 2007 
March 1.84 (UT dates are used throughout this paper; 
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iNakano & Itagakil [2007b . A spectrum obtained on 2007 
March 4.34 showed that SN 2007af was a SN la at least a 
week before peak brightness (Salg ado et al.ll200"7l) . The host 
galaxy of the supernova is NGC 5 584, an Scd galaxy wi th a 
recession velocity of 1638 km s" 1 ( Koribalsk ietaTll2004l) . 

2.1. High-Resolution Spectroscopy 

We ob served SN 2007af with the ARCES echelle spec- 
trograph (Wang et al. 2003) on the ARC 3.5 m telescope at 
Apache Point Observatory (APO) on 2007 March 10. We 
obtained two 1440 s exposures covering the spectral range 
3200-10,000 A at a spectral resolution of R w 33,000 and a 
signal-to-noise ratio (S/N) of 36 per pixel at the wavelength 
of the redshifted Na D lines. The data were reduced in IRAF y 
with the ECHELLE package using standard procedures. 

We also obser ved SN 200 7 af with the high-resolution spec- 
trograph (HRS; lTuUlll998t) on the Hobby-Eberly Telescope 
(HET) on 2007 March 31. The spectrograph was in its 
R = 30,000 mode, with a 2"-diameter fiber and the 316 £/mm 
grating centered at 6948 A, providing nearly complete wave- 
length coverage from 5095 to 8860 A. We obtained 4 spectra 
totaling 3100 s of exposure time and reached a combined S/N 
of 87 per pixel. The HRS data were reduced in IRAF with the 
echelle package. 

Final ly, we observe d SN 2007af with the HIRES spectro- 
graph dVogt et alj 19941) on the Keck I telescope on 2007 April 
7, over the range 3 150-6000 A. We obtained a single 900 s ex- 
posure through a 7."0 x 0"861 slit, yielding R w 48,000 and 
S/N = 47 per pixel. The HIRES data were processed with the 
MAKEE data reduction package. 

2.2. Imaging and Low-Resolution Spectroscopy 

SN 2007 af was the target of extensive photometric follow- 
up observations with the 0.76 m Ka t zman Automated Imag - 
ing Telescope (KAIT; iLi etaT1l2000t iFilippenko et al.ll2001l) . 
continuing for over 4 months since its discovery. No pre- 
explosion BVRI images are available for subtraction of the 
host-galaxy light, but the supernova occurred relatively far out 
in the disk, away from signific ant contamina ting features. We 
used the DAOPHOT package (IStetsonlll987h in IRAF to per- 
form point-spread function photometry of SN 2007 af relative 
to various field stars in the KAIT images, which were cali- 
brated on five photometric nights with KAIT and the Nickel 
1 m telescope at Lick Observatory. 

We also obtained low-resolution spectra of SN 2007 af with 
the Kast spectrograph (Miller & Stone 1993) on the Shane 3 m 
telescope at Lick on 2007 March 13, April 10, and June 14, 
which were reduced in IRAF following normal procedures. 

3. RESULTS 

3.1. Light Curve and Low-Resolution Spectra 

We display BVRI light curves of SN 2007 af in FigureQ] We 
fitted the photometric data with the latest ver sion of the mul- 
ticolor light-curve shape method (MLCS2k2; lJha et al.l l 2007l) 
to determine the parameters of the supernova. We find that 
the time of B-band maximum was 2007 March 14.76 (JD = 
2,454,174.26), with an uncertainty of 0.12 d. The derived 
line-of-sight extinction to the SN is Ay = 0.39 ± 0.06 mag, 

9 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 
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FIG. 1 .— BVRI light curves of SN 2007af, from KAIT data. Photometric 
uncertainties are indicated by the plotted error bars, which in most cases are 
smaller than the displayed data points. The MLCS fits are shown by the black 
curves. 
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FIG. 2. — Low-resolution spectrum of SN 2007af one day before maximum 
light (black curve), compared to a scaled spectrum of the prototypical SN la 
1981B at maximum light (red curve; data from Branch et al. 1983; telluric 
absorption is present at ~ 6830 A and ~ 7560 A). The cyan curve shows the 
spectrum of SN 2007 af four weeks later, for comparison. 

with an extinction law of R v = 2.98 ± 0.33 (the Milky Way 
foreground reddening is 0.039 mag; iSchleeel et all 1998). 
The distance modulus to SN 2007af is (32.06-5 logH /72) ± 
0.06 mag, giving the SN an absolute magnitude of My = 
(-19.28 + 5 logHo/72) ± 0.08. The luminosity/light-curve- 
shape parameter is A = -0.04 ± 0.02, and the MLCS2k2 re- 
duced x 2 value of 0.53 indicates an excellent fit. Photometri- 
cally, SN 2007 af appears normal in every respect. 

In Figure [2] we show the spectrum of SN 2007 af one day 
before maximum light, along with a template spectrum of SN 
198 IB at a similar epoch. Ana lysis with the Superfit spectral 
fitting code dHowell et al .1120051) indicates that SN 2007 af is a 
very typical SN la, comparab le to archetypal events such as 
SNe 1981B and 1989B (e.g., lFihppenkol[i997h . 

3.2. Sodium D Absorption Lines at High Resolution 

We display the high-resolution spectra of SN 2007 af around 
the host-galaxy Na D lines in Figure [3] At least two absorp- 
tion components are present. In the ARC and HET spectra, 
the shape of the blue wing of the absorption profile (as well 
as the residuals from two-component fits) suggest that a third 
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— ARC spectrum, t - -4.3 d 

— HET spectrum, t = +16.6 d . 
| — Keck spectrum, t = +23.7 d - 



5920 5922 5924 5926 5928 5930 
Observed Wavelength (A) 

FIG. 3. — High-resolution spectra of host-galaxy Na D absorption lines in 
SN 2007af. The HET and ARC spectra are offset by 0.6 and 1.2, respec- 
tively, in the vertical direction for clarity. The tickmarks labeled "A," "B," 
and "C" at the top of the plot indicate the wavelengths of the three absorption 
components, and the dotted gray lines show the wavelength range used to 
measure the EW for components A and C. Note that the Keck spectrum has 
significantly higher spectral resolution than the other spectra (see ij2}. 

component at slightly lower velocity might be blended with 
the stronger absorption line, and examination of the higher- 
resolution Keck spectrum confirms this. 

We fit the Na Di and D2 absorption lines in each spec- 
trum with three Gaussian components using a Levenberg- 
Marquardt least-squares fit algorithm. 10 Each absorption 
component is given its own depth, width, and wavelength, but 
we hold the wavelength separation between the Di and D2 
lines fixed at its known value and we force the D1/D2 depth 
ratio to be the same for each component. Including the contin- 
uum level, there are therefore 11 fitting parameters for 123- 
308 spectral channels (depending on the dispersion of each 
spectrum). We list the fit results in Table Q] 

While the equivalent widths (EWs) of components A and 
C appear to change somewhat across the three spectra, in 
fact these components are so badly blended in the ARC and 
HET spectra that the individual EWs are not very well de- 
termined. Instead of measuring the EWs of these compo- 
nents separately with the Gaussian profile fits, we directly 
integrate the spectra from 5921.50 A to 5922.22 A. The com- 
bined EWs of components A and C in the three epochs are 
139 ±6 mA, 137 ±3 mA, and 148 ±3 mA." We therefore 
conclude that there is no change in these absorption lines to a 
3cr limit of 18 mA over the course of the observations. We es- 
timate that the Na I column density for components A and C is 
9 x 1 1 1 cm" 2 (assuming a Doppler parameter of b w 8 km s" 1 ; 
see TableQ}, and the change in the column between 4 d before 
and 24 d after maximum light is no larger than 2 x 10 11 cm -2 . 
The strength of component B, which is reasonably well sep- 
arated from the other two, does not change within the uncer- 
tainties over the course of the observations, and we estimate 
a Na I column density of (3 ± 0.2) x 10 11 cm" 2 for this com- 
ponent. We place 5ct upper limits on the presence of addi- 
tional unresolved absorption components at other velocities 
of 15 mA, 8 mA, and 7 mA (9, 5, and 4 x 10 10 cm" 2 ) for 
the ARC, HET, and Keck spectra, respectively. In contrast, 
the time-variable features in SN 2006X reached Na I column 

10 We obtain equivalent results with Lorentzian or Voigt profiles. 
1 All Na I EWs in this paper are given for the D2 component. The D1/D2 
depth ratio (and therefore the D1/D2 EW ratio) is 0.64. 



densities of up to 10 12 cm" 2 dPatat et alj|2007al) . 

3.3. Calcium H&K Absorption Lines at High Resolution 

The APO and Keck spectra extend far enough to the blue 
that we also detect Ca II H&K absorption lines in both the 
Milky Way and NGC 5584. We fit the four sets of lines si- 
multaneously using the same technique as in 33.21 As with 
the Na lines, the host-galaxy absorption can be well fitted 
with two or possibly more components. Their velocities and 
relative strengths, however, differ significantly from those of 
the Na absorption components. The Ca absorption systems 
have heliocentric velocities of 1643 and 1613 km s , and the 
higher-velocity system is approximately 3 times as strong as 
the lower- velocity system (EWs of 160 ± 4 and 49 ± 3 mA in 
the Ca II K line, respectively). These results suggest that the 
various absorbing clouds along this line of sight have differ- 
ent abundances and/or ionization states. Although the higher- 
velocity component appears somewhat deeper and the lower- 
velocity component appears narrower in the Keck spectrum 
than in the APO data, we do not detect any statistically signif- 
icant changes in the absorption line equivalent widths between 
the two epochs. 

3.4. Limits on Ha Emission 

The ARC and HET spectra cover the expected wavelength 
of the redshifted Ha line. We detect an Ha emission line in 
both spectra with a velocity of 1641 km s" 1 and a full width 
at half-maximum (FWHM) of 16 km s" 1 . In the earlier, but 
lower S/N, ARC spectrum (f = -4.3 d) we measure an EW 
for the Ha emission of 45 ± 1 1 mA, and in the HET spectrum 
(t = 16.6 d) we measure an EW of 87±4mA. At these epochs, 
the supernova had /?-band magnitudes of 13.32 ±0.02 and 
13.85 ±0.02, respectively. Within the uncertainties, there- 
fore, the Ha flux did not change between the two observa- 
tions. While Ha emission could be a signature of CSM, we 
also find weak [N II] A6583 and [S II] A6717 features at the 
same velocity. The presence of these additional lines, the lack 
of variability in the Ha flux, and the close agreement between 
the host-galaxy velocity and the Ha velocity suggest instead 
that the emission is coming from a nearby H II region. 

Assuming that the detected emission lines come from the 
host galaxy rather than the supernova, we place an upper 
limit on the Ha emission from the SN of 22 mA (5cr limit 
for a line width of 50 km s" 1 ) at t = 16.6 d and 33 mA 
at t = -4.3 d. Although the high-resolution spectra are not 
flux-calibrated, we can estimate the Ha limits in flux units 
by scaling the continuum flux level to the observed /?-band 
magnitudes. We calculate an upper limit to the Ha flux of 
- 2 x 10" 16 erg cm" 2 s" 1 (Lu = 2 x 10 37 erg s" 1 ). Using 
the model of [Cumming et al. (1996]), the corresponding up- 
per limit on the mass-loss rate in the progenitor system is 
not very restrictive (< 10~ 4 5 M yr" 1 ). Several other nearby 
SNe la have comparable or better limits from Ha observation s 
(ICumming et al]ll996tlMattila et al.ll2005HPatat et alJl2007al) . 

4. DISCUSSION AND CONCLUSIONS 

Our observations of Na absorption lines in 33.2l can be used 
to estimate the amount of hydrogen along the line of sight to 
SN 2007 af. The total measured Na I column density is 1 .2 x 
10 12 cm" 2 . The corresponding hydrogen column is (8.1 x 
10 17 ) /A" cm" 2 for a solar N a abundance of 12 + log(Na/H)= 
6.17 dAsplund et al.ll2005l) . where X = N(Na I)/N(Na) is the 
Na ionization fraction. If this material were all located in 
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the SN progenitor system in a thin shell w ith a radius of 
~ 10 16 cm (as suggested by lPatat et al.ll2007al for SN 2006X), 
the shell mass would be (8.5 x 10~ 7 )/X M . 

However, because of the lack of variations in the column 
density of the detected absorption lines and their close agree- 
ment with the host-galaxy velocity, the absorbing gas is more 
likely associated with the interstellar medium of NGC 5584. 
In that case, the relevant calculation for the mass of the CSM 
is based on the upper limits for additional absorption com- 
ponents. Using the 5 a upper limits on the Na I column den- 
sity of 9, 5, and 4 x 10 10 cm" 2 for the three spectroscopic 
epochs and the assumptions given above, we find correspond- 
ing upper limits on the shell masses of 6A/X, 3.6/X, and 
2.9 /X x 10~ 8 M Q , respectively. For comparison. iPatat et al.l 
d2007al) estimated a shell mass of 7.1 x 10" 7 M for SN 
2006X with the same model. If the circumstellar Na ions 
have not yet fully recombined by the time of our observations, 
then the true CSM mass could be significantly higher, but for 
SN 2006X substantial recombination had occurred within 2 
weeks of maximum light. 

There are two primary ways to interpret the absence of vari- 
able absorption features in SN 2007af. First, the progenitor 
system may diff er from the red-gian t companion in a recurrent 
nova model thatlPatat et al.l d20 07a) proposed for SN 2006X. 
In particular. lWang et al.l (l2007l) found that SN 2006X has the 
highest expansion velocity ever measured for a SN la, which 
raises the possibility that varying Na I D absorption features 
(and hence the proposed recurrent nova progenitors) could be 
associated with the s ubgroup of high velocity gradient SNe la 
dBenetti et al.1 120 05). If this conjecture is correct, it would 
suggest that the companion to the SN 2007 af progenitor was 
either a main-sequence star, a subgiant, or another C/O white 
dwarf, leaving the progenitor system relatively free of circum- 
stellar gas and dust. 

Alternatively, the observational differences between the 



Na D lines in SNe 2007 af and 2006X could be a result of the 
system geometry. If the mass lost by recurrent novae is con- 
centrated in the orbit al plane of the system, as appears to b e 
the case for RS Oph (lO Brien et all 120061: iBode et al.ll2007l) . 
then lines of sight that do not pass near the orbital plane would 
miss most of the circumstellar material. Such a configuration 
would allow the progenitor systems of SNe 2006X and 2007 af 
to be the same, as long as the viewing geometry is different. 

Monitoring of the Na D lines in a larger sample of SNe la 
will clarify the geometrical constraints, including possible ad- 
ditional geometries for the CSM. It will also allow us to ex- 
amine whether the presence or absence of variable absorption 
is related to any other observed properties of the objects. 
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TABLE 1 
Na D Fit Results 



Component A Component B Component C 



Spectrum 


Velocity" 


FWHM 


EW 


Velocity 11 


FWHM 


EW 


Velocity" 


FWHM 


EW 


X 2 




[kms- 1 ] 


[kms" 1 ] 


[mA] 


[km s" 1 ] 


[kms" 1 ] 


[mA] 


[kms" 1 ] 


[kms" 1 ] 


[mA] 




ARC 


1632.2 ±0.5 


11.6±2.1 


80 ±25 


1650.4 ±0.6 


9.5 b 


43 ±8 


1625.0 ±7.6 


26.9 ±7.5 


62 ±37 


0.74 


HET 


1633.2±0.7 


11.3 ± 1.2 


87 ±18 


1649.1 ±0.4 


12.3 ±0.9 


54±4 


1622.5 ±3.4 


15.3 ±4.1 


45 ±16 


1.68 


Keck 


1632.6 ±0.2 


10.1 ±0.4 


104 ±4 


1649.4±0.1 


7.3 ±0.3 


47 ±3 


1621.7 ±0.4 


8.5 ±0.8 


34±4 


1.33 



" Velocities are in the heliocentric frame. b According to the fit this line is unresolved, so we cannot establish an uncertainty on the line width. 



